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A new methodology to determine kinetic 
parameters for one- and two-step chemical models 

By T. Mantel 1 , F. N. Egolfopouios 2 L. C. T. Bowman 3 


In this paper, a new methodology to determine kinetic parameters for simple 
chemical models and simple transport properties classically used in DNS of pre- 
inixed combustion is presented. First, a one- dimensional code is utilized to in- 
formed steady unstrained laminar methane-air flame in order to verify intrinsic 
features of laminar flames such as burning velocity and trmpe- -e and concen- 
tration profiles. Second, the flame response to steady and unsu, y strain in the 
opposed jet configuration is numerically investigated. It appears that for a well 
determined set of parameters, one- and two-step mechanisms reproduce the extinc- 
tion limit of a laminar flame submitted to a steady strain. Computations with the 
GRI-mech mechanism (177 reactions, 32 species) and multicomponent transport 
properties are used to validate these simplified models. A sensitivity analysis of the 
preferential diffusion of heat and reactants when the Lewis number is close to unity 
indicates that the response of the flame to an oscillating strain is very sensitive to 
this number. As an application of this methodology, the interaction between a two- 
dimensional vortex pair and a premixed laminar flame is performed by DNS using 
the one- and two-step mechanisms. Comparison with the experimental results of 
Samaniego rt al. (1994) shows a significant improvement in the description of the 
interaction when the two-step model is used. 


1. Introduction 

During the past ten years, direct numerical simulation (DNS) of turbulent react- 
ing flows has been widely utilized to obtain physical understanding and precious in 
formation for modeling purposes. The recent articles of Poinsot. (1996) and Poinsot 
e.t al. (1990) can be consulted for a review’ concerning DNS of turbulent reacting 
flows. Although any kind of model is needed to solve the Navier-Stokes equations 
for a non-reacting system, closures have to he provided in order to model transport 
properties of the different species and chemical reactions as well. These two aspects 
can rapidly lead to tremendous needs of storage rapacity and CPU time even for 
the combustion of simple hydrocarbons such as methane. As an example, the re- 
cent detailed mechanism proposed by the Gas Research Institute (GRI) for methane 
combustion requires 177 reactions of 32 species. This kind of chemical scheme ran 
only be used in the computations of one-dimensional problems such as the study 
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of strained laminar premixed flames (EgolfopouLs 1994a-b). We can, however, cite 
the two-dimensional numerical study of vortex- premixed laminar flame interactions 
performed by Hilka et al. (1994) using a detailed mechanism (17 species and 55 
reactions). Thus, in order to investigate turbulent flames propagating in the com- 
bustion regime of existing devices, the chemistry and transport properties have to 
be drastically simplified. Currently, a one- step irreversible chemical model is used 
to performed parametric studies of complex flows such as flame-vortex interactions 
(Poir.sot et al. 1991 ), three-dimensional decaying turbulence interacting with a pre- 
mixed flame (Trouve & Poinsot 1994), or a diffusion flame (Vervisch 1992). In order 
to take into account the highly diffusive behavior of some radicals, two-step mech- 
anisms have been used in numerical studies of turbulent diffusion flames (Vervisch 
1992) and flame- vortex interactions (Mantel 1994). The difficulty of these simple 
models is to find realistic transport properties and kinetic parameters which corre- 
spond to the studied medium. In the DNS code used by these various authors, the 
transport properties are modeled using a temperature dependence for the dynamic 
viscosity and constant Prandtl, Schmidt numbers, and calorific capacity. For the 
chemical models, the activation energy E a and frequency factor B for each reac- 
tion have to be estimated. Usually, a high activation energy in the range of 30 to 
60 kcal/mol is considered. For premixed systems, the kinetic parameter of these 
simple chemical models are chosen to match the laminar burning velocity Si alone. 
An a priori global art i vat ion energy must be taken high enough to be realistic b» 
low enough to reduce the number of grid points required to resolve the flame (ge. ■ 
ally the lowei limit of the range 30 to 60 kcal/mol). In fact, the asymptotic ana - 
of Williams (1985) shows that 6j ~ j~ l where tf is the thickness of the rca ion 
zone and i the Zel'dovich number defined by J — E a (Ts — T a ) / R° Tjf . Here R‘ 
and represent respectively the universal gas constant and the temperature of the 
fresh and burnt gases. Once the activation energy is imposed, the frequency factor 
is tuned to find the chosen laminar flame velocity. However, since an infinity of cou- 
ple ( B. E„) exists for a given value of Si, additional features of the laminar flame 
have to be verified. Thus, this technique has to be improved in order to predict 
other intrinsic characteristics of the flame such as concentration of reactants and 
temperature profiles, especially in the downstream end of the flame where reactions 
take place. 

The motivation of this study is to provide realistic kinetic parameters for one- and 
two step mechanisms classically used in DNS of premixed turbulent combustion. 
To do so, a new methodology allowing the determination of kinetic parameters 
is proposed. This methodology allows to verify the following quantities: (1) the 
laminar burning velocity, (2) the temperature and concentration of reactant (and 
intermediate species for the two-step mechanism) profiles, nd (3) the strain rate 
imposed to the flame in the opposed jet flame configuration leading to extinction . 

Points 1 and 2 are performed using the PREMIX code (Kee et al. 1994). which 
has been modified to accept artificial species, constant molecular weight, constant 
calorific capacities for all the species, and modified heat of formation to predict 
the adiabatic flame temperature. Point 3 is numerically investigated by studying 
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F»‘.;i BE 1. Omit erflow flame configuration 

the oounterflovr cvpjKwd jot dame configurati«>n (Egolfopoulos i994i. For these 3 
pwuis. computations using the GRI-roech 2.1 mechanism (Frroldarh et ei 1595) 
are perform** j and utilized as reference cases *5r comparison with one and t wo- step 
mechanisms. 

In order to validate this methodology or real configurations, the response of a 
premixed laminar flame to unsteady strain is numerically investigated in two differ- 
ent configurations using one- and two-step chemical models and simple trails j>or! 
properties: 

- the opposed jet flame submitted t an oscillating strain rate 

- tb*- vortex -premixed laminar flame interaction experimentally studied by Samanicgo 
rt ai (1996) 

In 'he first unsteady configuration, the effect of thermo- diffusive properties of the 
mixture is i fives* ated. It appears that this effect seems to have a strong influence 
on the unsteady behavior of the heat release rate. On tills configuration, both 
one and two-step models allow a good description erf the behavior of the flame, 
hi the case of the vortex -premixed laminar flame interaction, a significant effect 
of the diffusivitv of the intermediate species tin the heat release is observed when 
the- two-step mechanism is employed. Comparison with the experimental results of 
Samanicgo rt al. ( 1996) shows an improvement in the description of the interaction 
using the two-step model and simplified transport model. 

2. Lean premixed laminar flames submitted to a steady strain 

2. 1 Presentation of the. counterfiow flame configuration 

To study the ability of one- and two-step models to describe the resjumse of a 
laminar flame to stretch, the counterflow flame interaction configuration is chosen. 



152 


T- Mantel. F. M. Efolfopenlo.t & C T Bowman 



FlGl ir. 2- E'TiIntioiiof the laminar burning velocity and extinct »ou strain rate ver- 
sus the equivalence ratio :• : GRI-Oiech: : model: : E, = 35kral/niol: 

: E m - 60kcal/mol. 



FiGIRt: 3. Evolution of the activation energy in kacl/inol and frequency factor 

in cin.moi. s foe reaction ( 1 ) versus the equivalence ratio. : PREMIX code: 

: asympt;>tic analysis 

Such a configuration (see Fig. 1) has been widely studied both experimentally 
( Chung rt al. 1986. Law et al. 1986) rend numerically Fgolfopoulos 1994a.b). The 
main goal of these studies was to determine the extinction and flammability limits 
of laminar premixed flames Extinction strain rates and laminar flame velocity have 
be**n determined for a wide range of equivalence ratio for various air /fuel mixtures 
for premixed laminar. 

Here, this problem is treated using a odr solving die equatious of mass, inotueii 
turn, energy, and species along the stagnation streamline of the counterflow opposed 
jet flame configuration. Details concerning the equations and boundary conditions 
are given in Egolfopoulos (1994a). 

The conditions of our simulations arc those retained by Egolfopoulos (1994) 
The temperature of the unburnt mixture l methane-air ) is 300 K. and the distance 
separating tin* nozzles is 0.7 cm. 






Determination of kinetic pt.rameters for simple and complex mechanisms 153 


2.2 A one-step model for the combustion of lean methane -oxygen flames 

A new model for the combustion of lean methane-air flame is proposed following 
the methodology presented in the introduction. The global one-step reaction for 
lean methane-air combustion is: 


CH 4 + 20 2 - COj + 2H 2 0 (1) 

and the reaction rate for this reaction is expressed by 

RR = F[CH 4 )(02] 2 exp(-£./**T) (2) 

Computations are performed using the PREMIX code including reaction ( 1 ) and 
multi-component properties. The parameters B and E m are thus determined for 
each value of the equivalence ratio p. Figure (2) shows the evolution of St and 
K, wl for the global mechanism and for the solution obtained from the GRI-mech 
mechanism. The values for B and £. are presented in Table ( 1) as a function of o 
varying front 0.55 to 1. 

Taxt additional cases are presented in Fig. (2) by keeping constant kinetic param 
eters: (i) E m = 35 kevd/md; B = 4.1 10 21 (ii) E m — GQkcal/mol ; B = 6 10‘ 4 . These 
se:s (A parameters are determined for ^ = 1.0 and are kept constant for the other 
values of ♦. For these two cases, both St and h rr i are not correctly predicted. For 
some values of <*>, K tt , is even over-predicted by a factor of two (see Fig 2). 
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Tabic 1 . Kinetics parameters for global reaction defined by Eq. ( 1 ) used in Fig. ( 1 ) 

The evolution of £ as a function of <P obtained in the present study can be 
compared with the asymptotic analysis of Clavin ( 1985), who proposes an expression 
for the laminar burning velocity St as a function of B and 3. 
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where 


/ X m e’ x dX 

(4) 

§|7 


v V 'BC:-.C%exp(-3/a) 

(5) 


Here, v, J'„ W are the stoichiometric coefficient, the initial mass fraction, and the 
molecular mass of the deficient species. The molar concentration of the oxidant and 
the fuel in the fresh mixture are denoted by C"* and C*’ u . In Eq. ( 4 ), n is the 
order of the reaction and .Y a variable of integration defined by X = d(l - 0) where 
6 — {T - r.)/(T» - T m ) represents the reduced temperature, 
hi the case of the global reaction (1), Eqs. (4) and (5) become: 

T = 2 (6) 

— = » , CH 4 rr^ s -BCo^ „Cci# 4 , *cxp(-$/a) (7) 

T r 

To estimate the thermal diffu&ivity in the burnt gases, we use the classical relation: 



with b = 0.76. 

Reporting Eqs. (6-8) into Eq. (3). we obtain: 
t PrW 2 

B = I 5 ^V2 — -(l-o^-^expld/o) (9) 

Due to the assumptions used in the asymptotic analysis ( constant calorific capaci- 
ties. thermal, and species diffusivities ). Eq. (9) constitutes a first approximation Cor 
B. The values for B given by Eq (9) are compared with the result obtained using 
PREMIX. Asymptotic analysis exhibits higher values for B compared to PRE- 
MIX. This is also noticed by Rutland (1989). who studied the propagation of a 
one- dimensional premixed laminar flame using a one-step chemical model and by 
considering constant transport properties. 

t.S Kinetics parameters of one- & two-step models for Icon premixed laminar flame 

In this section, the kinetic parameters for one- and two-step models are deter- 
mined fix the combustion of a methane-air premixed laminar flame with an equiva- 
lence ratio of 0.55 using the PREMIX code. Particular conditions fix the transport 
properties are considered The dynamic viscosity is expressed according Eq. (8) 
and constant Prandt) and Schmidt numbers are assumed. Calorific capacity is also 
assumed constant and the molecular weights of all the sjwcies are equal. To do 
so, the PREMIX code had to be modified to accept artificial species and modified 
transoort properties. 

The motivations of these choices are directly related to the DNS code applied to 
complex flows such as vortex- premixed flame or turbul«*nce-preniixed flame interac- 
tions. 
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2.3.1 The one-step model chemical model 

In this model, the chemistry is described by a single step irreversible reaction: 

A{ react wits) — » P( products) (10) 

The reaction rate of this reaction is expressed using a classical Arrhenius law 

«M* BCaexpl-EJILT) (11) 

For this simplest chemical model, 4 parameters appear: B. (A If), and Le * 

((AW) being the heat released by the reaction). Since He* and (AW) can easily 
be determined (by using binary diffusion coefficient for the Lewis number and by 
matching the folly burnt gas temperature for (AW)), we have to determine B and 

E.. 

2.3-2 The two-step chemical model 

The two-step mechanism initially proposed by Zd'dovkh ( 1948) consists of a first 
order chain branching reaction and a second-order termination reaction: 

A + X -* 2X (12) 

X+X-P (13) 

The use of a two-step mechanism significantly increases the number of unknowns 
Now. S parameters have to be determined: B| , Bj, £«,* E mj . (AW)|. (AW)?. Ic a, 
and Lt\ where (AW)i and( Aff)j represent tbe heat r ele as ed by the first and by 
thc second reaction. To reduce tbe number of unknowns, some realistic assumptions 
can be proposed: 

- the first reaction has a high activation energy and is thermo-neutral (Linan 1974) 

- the second reaction has a zero activation energy and liberates all the heat ( Linan 
1974 ) coefficients 

These assumptions lead to simplified expression for the reaction rates of the 
reactions (12) and (13). 


WW, =W,C A C x exp(-£:. 1 /W.r) (14) 

Wf?2 = (15) 

Moreover, since the H atom provides a crucial source of radicals and plays a deter- 
mining role in the submechanism H 2 — O 2 (Glassman 1967), we relate the interme- 
diate species of the two-step mechanism to the H atom Thus, from binary diffusion 
coefficients, the Lewis numbers for A and X are: Lc* = 1.0, Lf \ — 0.15. From 
these considerations. 3 parameters still have to be determined: By . W 2 , and E„ { . 

To determine the remaining unknowns of the one- and two-step mechanisms, the 
methodology previously described is applied. 
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FKURE 4. Tcn:|Mr*t.iie and A uta^ fraction obtain'd with the one-step model 
{reduced by its value iu the fresh gas) at the trailing edge of the flame. o : GRl- 

tuech (CH» mass fraction): : E m = 30kcal/iu<»l B — 1.82«8: : E n = 

45k>'al;'m<>! 4 93< 1ft F, = COkral/tmil B = 1# i3 i B i'i mot. si. 



X |oa| x | cm) 

FlGl HE 5. R'duced aifthan* 1 aud H cone: nitration and temperature profiles aer«jss 
the laminar Haim* front. : GRI-niech: : two step uicrlianisin. 

2.3.3 Result.* 

First, the influence of kinetic parameters of one- aud two-step models on tlie flame 
structure is analyzed. The flame structure is very sensitive to the couple ( B. E„ ) 
especially in the trailing edge of the flame (s«-e Fig. 4). Since both B ami E a vary, 
it is difficult to know which of these two parameters influences the gradients of 
temperature and concentration. The profiles of $( «, = V< n,/lcn,.« in the burnt 
gas side seems to be very critical in the opposed jet r< Hibernation. When the flames 
interact between them, incomplete combustion by leakage of the fuel can lead to 
sudden extinction. This is jiarticularly true for the case E„ - 30 keal/mo!. B — 
1.82 10 s mt»l~ , s - ' for which the spreading of the CH| profile is more pronounced. 
For these values, the one step model predicts an extinction strain rate of 
whereas the experimental results give I\ r .-p — 200s 1 i Egolfojmulos 1904a). 

This ran he explained by noticing that for high activation energy, the thin reaction 
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FlOlRE 6 . Evolution of the heat release integrated across the flame front ( reduced 

by the unstrained value) in function of the strain. : GRI-rnech: : one- 

step mechanism; ; two-step mechanism. 



FtOi'RE 7. Evolution of the strain rate in function of time for different frequencies. 

; 1 Hz; ; 40 Hz; : SO Hz. 

zone is located at the downstream end of the temperature and concentration profiles 
(since the reaction zone is proportional to /J -1 ). Thus, for larger values of E». 
low strain rates only affect the preheat zone. As the strain rate increases, the 
temperature and concentration profiles are steeper, and the reaction zone starts to 
be affected by the strain. This effect is emphasized in the twin flame configuration 
where the distance separating the two reaction zones is a key parameter in the 
processes leading to extinction. The determination of the extinction strain rate 
(in ihe opposed jet configuration) depends directly on the good prediction of the 
position of the reaction zone in function of the inlet mass flow rate and, consequently, 
the strain rate. 

Figure ( 5 ) represents the flame structure using the two-step mechanism. Here, the 
concentration of the intermediate species is also of interes* because of the quadratic 
dependence on V‘x on the heat release rate (see Eq. 15). The maximum value of X 
is chosen by matching the maximum value of the H atom concentration given by the 
GRI mechanism. Once B\ and E, t are chosen to match the & and ♦ 4 profiles, the 





158 


T Mantel. F. N. Egolfopoula .» ft C T. B own an 


maximum of ♦ \ is directly related to the frequency factor of the second reaction 
B>. We also observe a sjiace shift of the coni|»ared r«» the 4</y profile This could 
he overcome by decreasing B\ in order to change the production of X. However, 
tile decrease of B\ directly leads to a thickening of the fiante in the burnt gas side 
of the flame and. consequently, to a different response of the flame to strain. Thus, 
a compromise between the profiles of 6. $ and has to be found in order 
to obtain the right extinction strain rate in the opposed jet flame configuration. 
Finally. Fig. (6) shows the flame response to a steady strain rate obtained for the 
one- and two step • .vdels and the comparison with the solution given by the GRI 
mechanism. Heie. the heat release rate integrated across the flame (normalized 
by the unstrained value) is presented Both the one and two-step mechanisms 
allow «s to find the correct extinction straiu rate (within 10'/t of error). We also 
notice that the Lewis number effect is also observed by using the simple transjrort 
properties described in section (2.3). Since the Lewis mnnl>cr based on the limiting 
species (here CH* ) is less than unity (Lf.\ = 0.95). a positive stretch applied to 
the flame increases the heat release As the stretch increase-., the reaction zones 
are pushed toward the stagnation plane, and reaction cannot lie sustained due to 
shorter residence time (Law 1988). 

3. Lean premixed laminar flames submitted to an oscillating strain 

y. 1 Analytic of the flame rrnpon.it 

The response of a laminar premixed methane-air flame to unsteady straiu is 
numerically studied using one and two step chemical models. The unsteadiness of 
the flow is obtained by imposing a sinusoidal velocity field at the inlet boundaries. 
The amplitude of the velocity variations is 20 '/i of the mean value inlet velocity. 
Three different frequencies for the velocity fluctuations are studied (1. 40. and 
SO Hz). The strain rate applied to the flame varies from 70 to 1G0 s -1 . corresponding 
to Karlovitz number varying from 0.33 to 0.75. Here, the Karlovitz number is 
defined by A’« = r r A where A’ is the strain rate anti r r a chemical time scale 
defined by r c — o,/Sj . o „ being the thermal diffusivity in the fresh gases. The 
order of magnitude of the Karlovitz numlier is typically representative of the fiamelet 
regime defined by the Klimov- Williams criteria (An < 1). 

The time evolution of the heat release rate integrated across the flame (non- 
dimensionalized by tin* unstrained value) is presented Fig. (8). As a first obser- 
vation. no phase shift is observed between the one- and two-step chemical models 
and the solution given by the GRI mechanism. The slight asymmetry between the 
slopes corresponding to the extension and relaxation observed by the GRI mecha- 
nism is also described by the two-step model. Figure (9) represents the iieat release 
amplitude (normalized by its steady strained value) for different frequencies. The 
amplitudes of the fluctuations are underestimated by the simple models even if the 
tendency is well reproduced. As previously observed by Egolfopoulos (1994a). at 
low frequencies the flame behaves like in the steady case whereas at higher frequen- 
cies. the amplitude of the fluctuations decrease. The attenuation of the heat release 
amplitude at higher frequencies is explained by the fact that the disturbances are 
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FIGURE 8. Evolution of the heat release integrated across the flame front {reduced 

by the unstrained value) in function erf time ( Le a — 0.95). : GRI-mech; 

: one-step mechanism; : two-step mechanism. 



i 


(Mai 


Figure 9. Evolution of the heat release amplitude (normalized by its unstrained 
value) in function of the frequency. © : GRI-mech; + ; one-step mechanism; o ; 
two-step mechanism. 
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rapidly attenuated by diffusion effects (Egolfopoulos 1994a). 

S.t Response of a flame wttk unity Lewis number to an oscillating strain 
We also notice that the heat release rate given by the two-step mechanism is 
always lower than the results given by the one-step mechanism. This can be ex- 
plained using the asymptotic analysis of Seshadri k Peters (1983) who studied the 
structure of a planar premixed laminar flame submitted to stretch. Considering 
a high activation energy for the first reaction, the authors derived an asymptotic 
expansion for the temperature. They found that the first order temperature can 
be expressed as a function of stretch and Lewis numbers for the reactant and the 
intermediate species: 



Le A ~ 1 

Ua 


+ 



l - Le\ 

Le x 


/o l 


+ 


2X1 

Le x I 


(16) 


The subscript 0 refers to the axial coordinate where Yx is maximum, / 0 is a 
function always positive, A‘ is a non- dimension alized stretch, and (—AHj) is the 
non-dimensional ized heat of reaction of the recombination step. The relation (16) 
points out the respective roles of the diifusivities of the reactant and of the interme- 
diate species. Considering only the first term on the RHS of Eq. (16), for positive 
stretch the temperature increases for Le A < 1 . For Le A = 1 , the temperature 
remains constant equal to the zero order temperature regardless the value of the 
stretch. This recovers the classical conclusions of the role played by the Lewis num- 
ber of the reactant on the dynamic of stretched flames (Clavin 1985, Law 1988). 
The second term on the RHS of (16) enhances the effects of diffusivity of the inter- 
mediate species on f he dynamic of stretched flames. Since radicals are mostly very 
light species, they have high diifusivities leading to Lewis numbers significantly less 
than unity (here Le \ 0.15). Thus, in the case of positive stretch, the diffusivity 

of the intermediate species tends to decrease the temperature and, consequently, 
the local laminar flame speed. This result points out that even for Lt x = 1 the 
flame can be sensitive to stretch effect and exhibits local variations of the laminar 
flame speed not only due to compression of the reaction zone. 

Moreover, under some circumstances, a positive stretch can produce a decrease 
of the heat release rate when the Lewis number of the reactant is slightly less than 
unity. This is observed in Fig. (8) for the frequency 1 Hz and Le A = 0.95 where 
the normalized heat release rate goes under unity. 

In order to characterize the effect of a slight variation of the Lewis number, the 
response of the flame to unsteady strain is analyzed by imposing the Lewis number 
for the reactant equal to unity. This slight variation of Lt\ has a strong consequence 
on the flame response. Figure 10 shows the evolution of the heat release integrated 
across the flame front for Le A = 1.0 and for 1, 40, and 80 Hz. The results issued 
from the one and two-step calculations are in opposition of phase compared to 
the solution given by the GRI mechanism. Due to the compression of the reaction 
zone, the integrated heat release is less than unity for both one- and two step- 
models. Moreover, for the two-step mechanism, the diffusion of the intermediate 
species also contributes to the decrease of the heat release as previously explained. 
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Figure 10. Evolution of the heat release integrated across the flame front (reduced 

by the unstrained value) in function of time ( Lc^ = 1.0). : GRI-mech; : 

one-step mechanism; : two-step mechanism. 

This behavior is well summarized on Fig. 11, in which the scatter plot of the heat 
release rate versus the strain rate is represented for all the frequencies ( 1 , 40, and 
80 Hz). Very clear correlations are observed, and different signs for the slopes are 
found between the cases Lca = 0.95 and Lt^ = 1.0. 

This seems to indicate that the thermo-diffusive properties of the mixture is a 
first order parameter in the behavior of strained laminar flames. 

4. Vortex-premixed laminar flame interaction 

The configuration investigated here concerns the interaction between a two-dimensional 
vortex pair generated by acoustic excitation and a V-shaped air- methane premixed 
laminar flame stabilized on a heated wire. A counter-rotating vortex pair propa- 
gating itself by mutual induction interacts with an initially planar premixed flame. 
Figure 12 shows the vorticity and heat release fields during the interaction ( t = 5ms). 

Here, the Lewis number based on the reactants is taken equal to unity. This prob- 
lem has been extensively studied both experimentally (Samaniego et al. 1996) and 
numerically (Mantel 1994). Here, a lean methane-air flame is investigated (equiva- 
lence ratio = 0.55). The initial conditions for the simulations are obtained from the 
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FIGURE 11. Scatter plot of the integrated heat release across the flaine front 
(reduced by the unstrained value) vs strain rate. 


experiment. The characteristics of the interaction are Vp/Si = 66.8, */df = 25.7. 
and D/Sj — 104.8, where I p, Si represent respectively the displacement velocity of 
the vortex pair and the laminar burning velocity: and s. A/, and D are the distance 
between the center of the vortices, the laminar flame thickness, and the distance 
separating the vortex pair from the laminar flame. 

Details concerning the geometry and liagnostic techniques can be found in Samaniego 
et al. ( 19°G). Information concerning the equations solved in the DNS code and 
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FIGURE 12. Vorticity (top) and heat release rate (bottom) fields at t = 5ms of 
the interaction (from Mantel et al. 1996). 

the computational configuration are presented in Mantel (1994). 

Since radiative heat looses effects have been found negligible during this interac- 
tion (Samaniego 1996, Mantel 1994), adiabatic conditions for the flame are taken 
for the simulations. Figure 13 shows the time evolution of the flame length (non- 
dimensionalized by its initial length). Comparison with the experimental results of 
Samaniego (1996) points out that the dynamic of the interaction is well reproduced 
by the simulations. The time evolution of the minimum heat release rate integrated 
along a normal to the flame and encountered along the flame is also shown Fig. 13. 
As long as the interaction goes on, the vortex pair increases the flame length and 
leads to a decrease of the heat release rate at a location in front of the vortex 
nan This is qualitatively well described both by one- and two-step mechanism. 
In this configuration, the two-step model allows a significant improvement in the 
dc; cription of the decrease of the heat release rate. 

5. Conclusions 

This paper presents a new methodology to determine kinetic parameters of one- 
and two-step chemical models classically used in DNS of premixed combustion. By 
using a one-dimensional code in which simple chemical models and simple transport 
properties are implemented, the kinetic parameters are determined in order to verify 
( 1 ) the laminar burning velocity, (2) the temperature and concentration profiles, and 
(3) the extinction strain rate of t, laminar flame in the opposed jet configuration 
(counterflow flames). To do so, the results issued from these simple models are 
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FIGURE 13. Time evoluion of the flame length (top) and minimum heat release 

rate integrated across the flame. : one-step mechanism; : two-step 

mechanism. 


compared in detail with results obt ..ned from the GRI-mech 2.1 mechanism and 
multi-component transport prope ties. 

Applications of these simple mechanisms and transport models on two unsteady 
configurations show good behavior of these models. In the case of a flame submitted 
to an unsteady strain, both one- and two-step models describe qualitatively wt>l the 
dynamic of the flame and the heat release amplitude for different frequencies. In 
this configuration, no obvious improvement is obtained with the two-step model. 
However, the Lewis number based on the reactant seems to be a determining pa- 
rameter for laminar strained flames. A slight variation of the Lewis number from 
0.95 to 1.0 leads to a completely different behavior of the flr . . 

The interaction between a two-dimensional vortex pair wi:h an initially planar 
premixed flame is also analyzed by DNS using the one- and two-step chemical 
models. In this case, comparisons with the experimental results of Samanicgo et ai 
(1996) shows that the two-step chemical model allows a better description of the 
interaction. However, further work is needed to investigate the effects of a slight 
variation of the Lewis number on the behavior of the flame during the vortex-flame 
interaction. 

This methodology can be improved by studying a configurat ion of laminar strained 
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flaute more representative of turbulent premixed flames. This concerns the conn 
terflow flame configuration, but with hot products ou one side and reactants on 
the other side Further work is in progress to « xamine tlie I avior of simple and 
transjvort model*. 
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